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I - LOW TEMPERATURE QUARTZ CRYSTAL OSCILLATOR 


GENERAL PRESENTATION 

The frequency stability of quartz oscillators is affected by 

- the noise of the electronics : thermal noise and 1/F noise 

- the peculiar frequency noise of the resonator : frequency random walk 
noise and 1/F noise. 

When the noise of the electronics is reduced at the lowest possible level 
compatible with the state of the art of electronic components, the limit of sta¬ 
bility of the oscillator corresponds to the resonator noise. 

- the frequency random walk is essentially due to temperature fluctua¬ 
tions and can be reduced by using high quality temperature controlled 

( 2 ) 

ovens 

- the 1/F noise level indeed defines the ultimate limit of stability. To 
the 1/F spectrum of the frequency domain corresponds a floor, the flicker floor, 
in the time domain. 

It was known that at very low temperature the internal losses and the 
acoustic attenuation of quartz material are reduced and therefore the Q-factor 
of resonators is increased. 

The present work shows that 

- the i/F noise level at room temperature is correlated with the Q-factor, 

4 

following a 1/Q -law. 

- at low temperatures it is possible to minimize, by at least one order 
of magnitude, the 1/F noise level. Improvement of the frequency stability 
follows. 
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The study of the noise levels as a function of temperature required 

- cryogenic equipment operating at 1°K 

- test set for crystal adapted to very high Q's and low motional resis¬ 
tances 

- frequency noise measurement system for pairs of resonators at 300°K am; 
1 °K. 


1/F noise is an universal noise, which can be observed in almost all phy¬ 
sical phenomena and systems. Its origin and mechanisms are not known. A theore- 

4 

tical study has been undertaken, and the first task was to explain the 1/0 de¬ 
pendence law, which was experimentally demonstrated. 


CRYOGENIC EQUIPMENT 

A liquid helium dewar with double isolation was used. Its capacity is 5 
and gives an autonomy of one or two days, depending on the number of leads and 
cables sunk in it. This small size was choosen in order not to expend to much 
liquid helium when the different quartz samples successively are cooled and 
heated. In the second phase of the project a 30 liters dewar will be used to 
maintain Xtal oscillators below 4K during periods of time of the order of seve¬ 
ral weeks. Vacuum pumps enable to reach temperatures between 4.2K and IK. Tempe 
rature is measured by means of carbon and germanium resistors, and also by mean; 
of the bath pressure when below 4.2K. 

For guide measurements, the resonators are directly sunk in the helium 
bath. If more stability, mainly versus temperature, is required they are enclo 
sed in a copper can. 

This cryogenic equipement was ready to use in February 1980. 
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CRYSTAL PARAMETERS AT l,OW TEMPERATURE 


The measurements were performed by using the * transmission method. 


Measurement system 

Measurements at low temperature present some difficulties because of the 
length of the cables between the resonator at the bottom ol' the dewar and the 
source and detector which are at room temperature. Therefore it is necessary 
to place the * network with the resonator in the bath. 

The test set is represented on fig. 1. 
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The n network must be adapted to 50 Cl , the characteristic impedance of 
the cobles. But. the vnlues of the parallel resistors must be chooses as ,i (uiw- 
tion of the motional resistance R^ of the resonator. This last one can become 
close to 1 Cl at very low temperature. Therefore two » networks were used. Theii 
values are given on fig. 2a and 2b. 



The n network head must be constructed with care. To avoid spurious con¬ 
tact resistances at temperatures below 4K all the electrical contacts were gold 
platted. The n head is represented on fig. 3. 


pins to 
connect 
Xtal 



Frequency, motional resistance and Q-factor were measured between IK and 
20K. This temperature range corresponds to the diminution of the acoustic atte¬ 
nuation, the phonon-phonon interaction peak being located at 20K, 

The highest Q-factor was obtained with G4 and corresponds to Q = 150.10^ 
at T = 1.5K. 
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and 5. 


The Q and R values of the resonators provided by RADC are shown on fig. 4 

The corresponding frequency shifts are given on fig. 6. They are of the 
order of 5 to 6 KHz. the temperature coefficient of the frequency was also meac 
red. It decreases with temperature, as expected (fig. 7). 

At room temperature, the O-factor values of these resonators ar<> clos'd, 
hut at low temperature a large dispersion is observed. The inverse O-factor can 
be written 

1/0 = 1 / 0 . 1/0 + 1/0 
b s e 

This relation points out the different loss contributions due to the crystal 

bulk (1/0, ), the surface (1/0 ) and the environment (1/0 ) constituted of the 
b s e 

electrodes and the mounting. 1/Q^ is representative of the crystal properties 
and depends on the phonon interaction processes. Its contribution is predominan' 
at 20°K and certainly also at room temperature. 1/Q^ decreasing with temperatur' 
as T 4 it becomes negligible at low temperature. The influence of 1/Q^ and 1/Q^ 
appears and leads to a threshold. These two last contributions are mainly repre 
sentative of the resonator technology. 



The 1/f frequency noise of resonators cannot be measured in an oscillatin/ 
loop, because the phase fluctuations of the loop due to the electronics are gent 
rally higher. The method which is used consists in driving in a it transmission 
network two identical resonators by means of a frequency source with high spec¬ 
tral purity (3)(4). j<j en tical resonators means resonators with identical frequer 
cies and 0-factors, or made identical with tuning capacitors and resistors. The¬ 
refore a preselection of pairs of resonators at the same frequency is necessary. 
The measurement system is shown on fig. 8. 
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When the two arms of the phase bridge are correctly balanced SO dB to 60 
iJB of rejection of the source frequency fluctuations are obtained. 

The resolution depends on the resonators 0-factors following the relation 
between the phase and frequency fluctuations d+ = - 2Q df/f , for Fourier fre¬ 
quencies within the resonator's bandwidth. 

On fig. 9 is shown a typical frequency noise spectrum. Three different 
noises appear : frequency random walk (1/F 1 ) for the lowest Fourier frequencies 
flicker noise (1/F) for the intermediate ones and white noise (F°) for the 

highest ones. 1/F 1 noise is related to the crystal temperature fluctuation AT 

^ ( 3 ) 

through its dynamic thermal coefficient a following the law Af/f = a AT/At 

This noise can be minimized by using thermostats with high accuracy. For AT cut 

\ _13 

a is 10 s/°C. This means that the achievement of a 10 frequency stability 

—8 

needs a temperature control as good as 10 °C/s. 

The white noise which is observed has nothing to do with the resonators, 
but is due to the measurement system itself. 


The level of 1/F noise was measured for a large number of pairs of resonn 
tors .it different frequencies, and with a large dispersion of the 0-factor va¬ 
lues. On fig. 10 are plotted the different noise levels measured on all those 
resonators and given at 1 Hz from the carrier. 


From these experimental data the dependence between 1/F noise and Q-facto- 
is determined by using a regression method (linear regression with the loganth 
mic scales). This leads to the general relation 


The values obtained for n and 1 depend on the selection of the experimen¬ 
tal points among the 18 available values. After mathematical treatment it can 
be shown that 3.8 4 n ■? 4.3. Therefore it is proposed to consider that the most 
probable value will be n = 4 . A small variation in the n value induces a drasti 
change of the constant X, and for the previous range of n the corresponding va¬ 
lues of x are 0.1 < X < 60. A reasonable fitting will correspond to x close 

to unity when n is made equal to the integer 4. 
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2) Low temperatures 


At room temperature there is no way for modifying the Q-factor of one re¬ 
sonator. An additional resistor will change the loaded 0, but not the unloaded 
one which is relevant for the noise measurements. The only possibility was to 
use different resonators with different unloaded Q-factors. 

But the acoustic attenuation is strongly dependent of temperature at low 
t empe-'ature;: and ran be either increased near the phonon-phonon interaction pea 
and near the impurity relaxation peaks, or decreased for temperatures below 10K 
as it has been shown. 

The previous noise measurement system, which requires a pair of resonat'- 
at the same frequency, is not aaapted for low temperature measurements. The 
frequency variations from 300K to 4K are of the order of a few KHz. But because 
of small differences in the crystal orientation, this frequency variation vers- 
sus temperature is not reproductible enough. Even when identical at room tempe¬ 
rature, the frequencies of a pair of resonator present differences of several 
100 Hz which cannot be cancelled out by using tuning external capacitors. 

Therefore a new measurement system for low temperatures was studied. Its 
schematic diagram is represented on fig. 11. 


F*f ISO 



(•'!«. 11 

1-13 





Each of the two resonators is excited by one of the side bands of an ampl 
tude-modulated signal (without carrier). At the output of the first mixer the 
intermediate frequency is affected by the total frequency noise of the resona¬ 
tors. The influence of the driving source is cancelled out by balancing the cha 
racteristics of the resonators. It is necessary to use a good frequency source 
for the modulating signal, even if its frequency spectrum in purified by divi¬ 
ding by 4000. 

The first mixer is followed by a low noise amplifier, with 20 dB of gain. 

A second mixing with a 2f reference frequency delivered from a digital 
phase shifter coupled to the divider gives a noise voltage proportional to the 
resonator's frequency fluctuations. A filter-amplifier cancells out the 4f fre¬ 
quency component. A FFT spectrum analyzer gives directly the S^(F) spectrum of 
the two resonators. 

One other advantage of such a measurement system, when compared to the 
previous one, is to enable easier and faster calibration, just changing by a 
small amount the modulation frequency and measuring the DC voltage shift at the 
mixer output. 

The 1/F noise of a pair of identical AT-cut quartz resonators (manufactu¬ 
red by CEPE) were measured at 300K, A.2K and at IK. The frequency difference 
between the two resonators waB 1 KHz. The noise levels for one resonator are 
given on table I. 


Temperature 

300 °K 

4,2°K 

1°K 

Q 

2.4.10 6 

5.10 6 

9.10 6 

S (1 Hz) 
y 

1.6.10 -24 

1.4.10 -25 

L - 

3.5.10 -26 


Table I 


These results confirm the dependence of 1/F noise on Q-factor and in this 

4 

case with the same resonator. The data does not follow strictly the 1/0 law. 

A comment on this point will be given in the next section. 
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theoretical approach 


It is known that the acoustic attenuation results from the interaction of 
the sound wave with the thermal phonons of the crystal. This is a fundamental 
phenomenon, which would appears even in a perfect crystal, i.e. without defects, 
because of the anharmonicities of the lattice, due to the nonlinear nature of 
the interatomic bonds. Also because of these nonlinearities are modifications 
of the phonon frequencies and of the sound wave velocity. 

Calculations of attenuation were made by using two different theoretical 

( 5 ) 

approaches. In the Landau and Rumer theory (1937) the sound wave is conside¬ 
red as a beam of phonons which are scattered by collisions with the thermal pho¬ 
nons. This theory does not take into account the interactions between thermal 
phonons and is valid only in the low temperature range, when the phonon life 
time and mean free path are large and therefore energy and momentum are well 
defined. At higher temperatures the life time becomes shorter. The uncertainties 
of energy and momentum increase and the selection rules arising from energy and 

momentum conservation laws break down. In the second theory due to Akhieser 

( 6 ) 

(1939) the sound wave is represented macroscopically (and not microscopi¬ 
cally as in the Landau-Rumer theory). The thermal phonon system is disturbed 
from equilibrium by the sound wave strain but tends to return to equilibrium 
because of the collisions between thermal phonons. In this case the phonons are 
localized and their mean free path has to be small when compared to the sound 
wave length. 


1) Harmonic approximation 

The main results of the lattice dynamics in the harmonic approximation 
will be briefly presented, (for more details cf. ref. 7). 


The potential energy $ 2 is assumed to be a quadratic function of the ato¬ 
mic displacements u ion < where n refers to nth Wigner-Seitz cell, a to the ath 
atom in the unit cell and i distinguishes the three components of the displace¬ 
ment vector 


l 

‘ ian 


i'a'n' 

6 

v ian 


u ian u l'o'n' 


(1) 


i * o * n ' 

are the second-order interactomic coupling parameters. 

ian 


1-15 



The kinetic energy of the crystal is 


V = (*> [ (2) 

ian d 

P. is the ith component of the momentum of the ath atom in cell n, and M is 
ion ° 

the mass of this atom. 

The Hamiltonian of the crystal is then 

h 2 = + * (3> 


The displacement u. and momentum P. can be written as a sum of waves 
ian ian 

with wave vector k, polarization j, travelling in the x(l) direction. 

Vi 2nik.x(n) 

U ion = J. [ 2NM> T k J) ] e ia (kJ) I e (4) 

. „ .. % 2» ik.x(n) 

P. an = -i [ [ hMa ^ ^i] e ia {k j) [ a(kj)+a + (-kj)] e (5) 

kj 


is the polarization vector, N the total number of unit cells. a + (kj) and 
o(kj) are the creation and annihilation operators. 


It can be shown that the Hamiltonian takes the simpler form 

H = l ti u ( k j) [ a + ( k j ) a ( k j ) + % ] (6) 

kj 

The sum over k is extended over the first Brillouin zone. 


The eigenvalues of the Hamiltonian H 2 gives the energy of the crystal 

E = [I N(kJ) ♦ '4] n u(kj) (7) 

kj 

which can be interpreted as the sum of the partial energies 

E(kj )={ N(kj )+/$]tiw(kJ) of each mode k and polarization j. N(kJ) is an integer. L 
teem of phonons it indicates the number of phonons'with energy "h*) (kj) in each 
partial wave. 
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2 ) Anharmonicity 


tn the harmonic approximation the various lattice waves, or phonons, are 
independant and there is no transfert of energy between them. Therefore a sound 
wave, which corresponds to a given lattice wave with a large number of phonons, 
will not be attenuated. 

a) Landau-Rumer theory 


Anharmonicities can be introduced by means of third (and/or fourth) order 
interatomic coupling parameters. In quantum mechanical model the anharmonic con 
tribution to the Hamiltonian takes the form 


H 


3 


3/2 

—H—. y _ 1 (-kj,k'j',k"j») 

3/p y L - -—--- 

2 ' 6N 7 kk'k" A>(kj) wfk"^') 

j j ' j" 


a + (kj) a(k'j 1 ) a(k"j") 


+ all other possible combinations of the a and a + operators. 


(a) 


These three phonon interactions must verify the energy and momentum consei 
vation laws. 


The sound wave can be considered as a lattice wave of mode KJ in which 
phonons are added so that their number N(kj) ia in excess of the thermal cquili 
brium number n(kj). The rate at which phonons are scattered out of thi3 mode 
enables to define a lifetime t(KJ) given by 


1 _ _1_ dN(KJ) 

1 (KJ) " “ N(KJ) - n(KJ) dt 


(9) 


If s(KJ) is the phase velocity of the wave the attenuation a is a = \/z is. 


The rate of change dN(KJ)/dt can be calculated by using the Golden Rule 
of the perturbation theory. 


The anharmonic coupling also perturbs the energy levels of the sound wave. 
This perturbation is calculated by using the standard second-order-time-lndcpen- 
deiTt perturbation theory. 

The validity of this method can be evaluated by considering thnt the the’ - 

mal phonons have a finite life time. Let t be lte average. In account of the 

Uncertainty Principle this finite lifetime makes the phonon energy uncertain 
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by an amount of the order of h/t. Therefore the method will not be valid as soor 
as the energy uncertainty is comparable to the sound wave phonon energy lift. Thu 
the Landau-Rumer method is valid when ffx « lift, or {It >> 1. This condition is 
satisfied for sound waves at high frequencies or at low temperatures. 


b) Akhieser Theory 

At higher temperatures the thermal phonon lifetime decreases. As a conse¬ 
quence their mean free path becomes smaller than the sound wave length X. The 
phonons can be considered as localized wave packets. The thermal phonon avera¬ 
ged frequency at temperature T is k T/1S (k_ : Boltzmann constant). If l is the 

D B 

wave packet length, the condition 1 <<X is equivalent to "tin << k n T. 

D 

In this case the method proposed by Akhieser consists in using a macros¬ 
copic model for the sound wave and the microscopic one for the thermal phonons. 
The anharmoniclties are introduced by means of a strain dependance of the phonoi 
frequencies. A detailed description of the method can be found in ref (8). 


The total energy of the crystal of volume V is the sum of the energy of 
the thermal phonons 

E = l (N(kj) + !0tl w (kj) (10) 

kj 

and of the mechanical potential energy due to the strains induced by the sound 
wave 

E .=C„Ti-+JiC„,n„n. + ... (11) 

sound aB aB aByfi aB oy 

C „ and C „ , have the dimension of elastic constants. 

OB GtpYti 

The derivative of the total energy with respect to strain gives the stress 


T = (a/a 

aB 


n ) 

aB 


th 


+ E .) 
sound 


(12 


" When in equilibrium the number of phonons in*mode kj is given by the Bose- 
Einstein distribution 


n(kj) =» (exp [Tj w(kj) / kgT] - l)" 1 


(13) 
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and this equilibrium corresponds also to T = n „ = 0. 

aB aB 


Inis enables to determine completely the stress 


T aB = C aBy* \« + J I t N < k J> " 


kj 


kj aB y« 


an 


aB 


(Id) 


The mechanical equilibrium equation 


B 

—r 

a t 


a t 


aB 


ax 


and the strain-displacement relation 


au 


a u 


n - ♦ —— 

aO ? ' 3x 3x 
B a 


(15) 


(16) 


enables to calculate plane waves solutions with wave vector K and frequency .fl¬ 
it 2*K.x-nt) 


u„ = u e (KJ) e 

B o B 

(17) 

i(2»K.x-nt) 

N(kj) = n(kj) + iN(kj) e 

(IS) 


In the absence of phonons the solution has the simpler form 

, ( i(2*K.x-Sl t) 

u = u e (KJ) e 0 

B o 6 


(19) 


The correction due to the phonons can be calculated by using a perturba¬ 
tion method. This leads to 


” ■ "o * |5r V« J tkjj 

o kj aB y6 

_ igL e (KJ) K i mm 

pVf! o ° 8 kj 3n aS U o 


( 20 ) 


The attenuation a is given by 

n 


a = - Im 


3 (KJ) 
o 


( 21 ) 
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and the new velocity of the sound wave 



These quantities are completely determined when AN(kj) is known. This cor 
responds to the calculation of the response of the phonon system to the sound 
wave, which is given by the Boltzmann equation 


3N(kj) 1 3o.(k,j) 3n ae 3N( k j) 3N(k,j) , 3N(kj) , 

31 " 2* 3n „ 3x 3k _ V 3x + 1 3t J coll 

aB Y Y a 


Remark : The term 
- 3t coll 

the collisions between phonons. 


is most important because it takes into account 
Without collisions and with tlr >> 1 and 


n << k^T/fl the present theory goes back to the Landau-Rumer theory. 


Three different types of collisions can be distinguished : 

- Normal processes (N), in which quasimomentum and energy of thermal pho¬ 


nons are conserved in three phonon collisions 

hk = hk' + hk" (24) 

flu) = ttu' + flui" ( 25 ) 

Total energy and total quasimomentum also are conserved. Therefore 

l ■»> 

l (3N(kj )/at 1 hk = 0 (27) 

kj C0X1 

The result of N processes is a distribution law which has the form 

[ftui(kj) - -nvk]/k T 

N N (kj) = (e 0 -U -1 (28) 


where 7 and T are such that the conservation equations are satisfied. 

- Umklapp processes (U) 

The quasimomentum is not conserved 

hk = hk' +■ hk" + hg (29) 

•flu> = flui* + flw" 
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(30) 




In the case of U processes, of N processes with some E and U processes 
(necessary to give a redistribution of the quasimomentum), or of E processes 
with some N and U processes (also to bring the phonon distribution to equili¬ 
brium) attenuation and velocity take the general form 



(35) 

(36) 


where y is an effective Grueneisen constant depending on the nature of the pro¬ 
cess involved, C the specific heat ; and t the mean relaxation time of the ther¬ 
mal phonons. 
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Introducing the low frequency young modulus E q = os^, the quantity AE=CTt j 

and. the relation between attenuation and 0-factor a = fl /(2Qs ), and between 

o o 

velocity and frequency (s-3 )/s = (tl-fi )/Si one obtains 

oo oo 


1 

0 


AE 

E 

o 


n 


o 


T 


1+n’T 1 

o 


(37) 


(9) 

In the case of a resonator fluctuation 4r in the relaxation time t will pro¬ 
duce a fractional frequency fluctuation y * 6(l/ft o the resonator frequency 
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(39) 


In terms of power spectral densities S^(F) and 
(37) in (39) one obtains 


(F) and using the relation 


.4 'aE' 


5 «t/t 


(40) 


4 

This relation shows the 1/0 law which was observed experimentally, bet¬ 
ween the noise levels and the Q-factors. 

For higher values of fJ Q T the Landau-Rumer method must be used. At the pre¬ 
sent state of this study the frequency change at low temperature was not calcu¬ 
lated yet. Only attenuation was evaluated. 

Near room temperature attenuation is proportional to T as shown by rela- 

4 

tion (35). At low temperature it becomes proportional to T . This corresponds 
to what is observed experimentally, even if the attenuation behavior is someti- 

4 6 

mes closer to T or T . The difference can be explained by the temperature depen 
dance of t, which does not appears explicitely in the relations. 

However temperature also appears in relation (40) in both 0 and E. When 

T decreases, 0 increases and therefore there is some compensation in the product 
4 2 

0 E . This can explain why the three experimental points of table I correspon¬ 
ding to 1/F noise levels at 300°K, 4.2K and IK does not strictly follow the 1/0 
law. 
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CONCLUSION 


The experimental results obtained during this study shown for the first 
time a correlation between 1/F noise and O-factor in quartz crystal resonators. 

This dependence was confirmed on the one hand by low temperature measure¬ 
ments and on the other hand by a theoretical model. 

This theory based on phonon interactions is a first approach to a more 
general theory of flicker noise. For this purpose it will be necessary to pre¬ 
cise the fluctuation spectrum of the phonon relaxation time, and calculate the 
frequency shift at low temperature by using the Landau-Rumer theory, in order 
to introduce completely the temperature dependence in the model. 

For a practical point of view, the large increase of Q-factor and the lo¬ 
wer noise levels open the possibility of achieving cryogenic quartz crystal os¬ 
cillators with better stability. 
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II - FAST WARM-UP SAW OSCILLATOR 


GENERAL PRESENTATION 


The characteristics of surface acoustic waves are modified by temperature 
variations. Those external effects are at the origin of the main causes of fre¬ 
quency instabilities in HAW oscillators. 

When the temperature fluctuations are slow enough, the thermal distribu¬ 
tion in the volume of the crystal is uniform. The usual model using a polynomia 
expansion up to the third order versus temperature is satisfactory to describe 
the static thermal behavior. For faster temperature variations a spatial distri 
bution of temperature within the crystal arises. Consequently, thermal effects 
involve temperature gradients which induce thermal stresses and strains. This 
is the "Dynamic thermal behavior"^ ^^^. These thermal stresses and strains arc 
coupled with the surface wave by the intrinsic and induced nonlinearities of 
the substrate. Velocity and frequency shift follow. 


The present work shows that 


- the fast temperature variations induce frequency changes which are pro¬ 
portional to the time derivative of the internal temperature. 


- the measured ST-cut dynamic temperature coefficient a is equal to 

-5 ST 

2.10 s/°K and is comparable to the coefficient of the AT-cut quartz crystal 

—5 

bulk wave resonator (1.3.10 s/°K). 


i ' 

- the measured dynamic TC's of doubly rotated cuts are larger in SAW than 


in BAW oscillators. 


The theoretical study was composed of 

- calculation of the temperature distributions and corresponding thermal 
stresses and strains by taking into account the mechanical mounting of the SAW 
device. 

<V 

- calculation of the dynamic T.C. a by using a perturbation method applle. 
to the general model of an acoustic wave propagating in a strained medium. 
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The measurements of the dynamic thermal behavicr required a temperature 
controlled oven operated by means of a computer. Numerical programs enable to 
perform given temperature cyclings, from the corresponding frequency responses 

M 

the a-coefficient are deduced. 


TEMPERATURE DISTRIBUTION 


A unidimensional model has been used. It consists in a crystal plate of 
2h-thickness and of infinite lateral sizes along a^ and (fig. 1). 



*(t) 


Fig. 1 : Temperature distribution 


A more accurate model should take into account the influence of the ford 
at the fixation points under the crystal thermal expansion. Only the bottom of 
the plate was submitted to temperature variations #(t). 

Let Gfa^.t) be the temperature inside the crystal. The boundary condition 

are 

0(-2h,t) = *(t) and 6(0,t) = 0 (1) 

The heat transfer with the surrounding medium on the propagation surface 
(» 2 = 0) ia considered to be very fast. 

Consequently, temperature diatribution within the crystal follows the gen 

(3) 

ral relation given by H. Carslaw 
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with 


/ “ p 
o ' 


2 2 , 
n ir kX 


-> (-1)“ $(X)dX 


where K is the thermal diffusivity constant of the crystal in the direction. 

Considering that «(t) is a slowly varying function of time, the integra¬ 
tion of I(t) has been performed by parts, and the time derivatives higher than 
the first order neglected. 


Then, the temperature distribution is given by 


a 1 / a \ 

O(a t) - - -1 *(t) + — ( - — + a h ) &(t) 

2h 3< \ 4h 1 ) 


$(t) is the first time derivative of ♦(t) 


*(t) = ♦(t),t 


The internal temperature 0(a_,t) having a spatial distribution, thermal 

* (4) 

stresses and strains appear. By following Holland , the adiabatic deformatio. 
n^j which takes place, is given by the following relations written in the matri 
cial reference system. 

(^]= A N a 2 + B N (6) 


with N s 1, 3, 5 


[n z ) = e(a 2 ,t) („’] - [C ZZ )- X [C ZN )[n N ] 


with Z = 2, 4, 6 


The different coefficients are given by the relations 


*» ■ li? '° h * e » 
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B-— (6 ~ + ~) 
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o = / 9(a ,t) da 
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0 

6=1 0(a 2 ,t) a 2 da 2 

-2h 


( 10 ) 

( 11 ) 

( 12 ) 


where and are the thermal dilatation constants and [C] are the elastic 
constants. 


By using the stress-strain relations 
mations and the temperature distribution, 
reference system) 


, tensions t.. 

ij 

are calculated 


induced by the defor- 
(also in the matricial 


' lt N ' ~ l ^ C NN ' “ ‘ C NZ^ C ZZ 


,-l 


t c z N l> 


<( V - t°N 


0(a 2 ,t) ) 


(13) 


and [t^l = ° (14) 

in account of the boundary conditions corresponding to a plate whic 
is free to expand. 


FREQUENCY SHIFTS 


Considering the temperature variations as a bias applied to the crystal, 
propagation properties are modified and in a natqral state coordinate system 
the nonlinear propagation equation is written 


p o u i,tt 


(A ikJm U j,m ) ,k 


(15) 


where p is 

o 

vibrations. 


the specific mass and u 


J 


the displacement due to the high frequency 


The boundary conditions corresponding to a stress free surface are 


A., . 

ikjm 


u . 

Jim 


0 


= 0 for 


(16) 



where A., . are the modified elastic constants, which can be written in the for 
ikjm 


A = C.. . +■)!,,. 

ikjm ikjm ikjm 


(17) 


where the H.. . tensor is considered as a small term with respect to the se 
ikjm 

cond order elastic constant C., . . It can be related to the temperature and the 

ikjm 

thermal stresses and strains in the following way 


H . = 6 t, + U. C, . + U C, . + h C. + a„ C.. . 0( a , t) (1R 

ikjm ij km i,p kpmj j,q kimq uv kmpquv C ikjm lkjm 2 


where t, , n and V, are respectively the thermodynamic tensions, the defor- 
km uv i,p 

mations and the displacement gradients. Here, a^ is the first order tempera 

ikjm 

ture coefficient of elastic constants for the considered-cut. 

In the case of SAW oscillators, the nonlinear coupling of the crystal witl 

the high frequency wave depends on its depth penetration. Then perturbation 

terms H~ are function of the step variable a„ and can be written following 
ikjm 2 

the relation : 


-- r (n) n 

ikjm L ikjm 2 


(19) 


(5) 

By using a perturbation method , relative frequency shifts are calcu¬ 


lated 
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Since the temperature distribution is proportional partly to *(t) and 
partly to the time derivative *(t), the perturbation terms H. k j m have two parts 


M • a #(t) + b.. , 4>(t) 

lkjm ikjm ikjm 


( 21 ) 
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Therefore, from relation (21) and the perturbation relation (20), the re - 
t.itivr frequency shift becomes after integration over the active surface of tie 


plate, 


—- - a$(t) + ai>(t) ( 22 ) 

o 


where the a^-term is related to the first order temperature coefficient of the 
frequency (static thermal behavior) and a is the dynamic temperature coefficiei: 
Theoretical results are presented in table I for different quartz-cuts. 


Table I 

Values of the dynamic thermal coefficient a for different orientations 

of quartz-crystal 


Crystal Cuts 

ST.X 

AT ,X 

Y.X 

((> - 6°20' 

FST 0 - -41“30' 

- 26° 

a in 10~ 6 S/°K 

0.04 

-1.4 

-3.5 

2.8 


1 

-35.0 

-87.5 

•>*4 

o 

o 

1 














TEMPERATURE CONTROLLED OVEN 


SAW delay lines built with different cuts of quartz crystal were submit¬ 
ted to temperature cyclings by using a temperature controlled oven . As show 
in figure 2, SAW devices and a LC-cut quartz temperature probe are mounted in 
the same oven. Temperature is measured with a quartz thermometer and controlled 
by a computer which drives the oven by comparing the measured temperature with 
a reference temperature. The temperature follows a variation law given by a pn 
gram included in the computer. Sinusoidal cyclings can be realized with an acc ■ 
racy of a few 0.001°C, as the SAW oscillator frequency is measured. 



Figure 2. Temperature controlled oven. 

D: SAW delay line; H: Heater; 
TP: Temperature probe; 

A: Amplifier; C: Directional 
coupler; DC; D.C. power supply; 
TH: Thermometer; FC: Frequency 
counter; COMP: Computer; 

PR: Printer; PL: Plotter; 

FR Freezer. 









KXPKHIMKfJTAL BKS'dl.TS 

Figure 3 shows an experimental curve obtained by sweeping temperature ov 
a temperature range of 4°C about 15°C with sweep frequency of 2.3.10 Hz for 

a reguiar ST-cut fixed in two points by epoxy on a flat heater. The measured 

_5 

dynamic temperature coefficient is equal to +2.10 S/°K and is significantly 

larger than the theoretical one. The oscillator frequency was 87.8 MHz. Ine sar 

plate was submitted to temperature cyclings applied on one end by a second heal 

source. Thus more inhomogeneous temperature distributions are introduced and 

the measured dynamic temperature coefficient obtained is ten times greater thar 

-4 

the previous one (-2.10 S/°K). 



13 14 15 16 17 

Temperature (®C) 


Fig. 3 -Experimental sinusoidal temperature cycling. Sweep 
frequency: 2.310 -<4 Hz, ST-Cut AO: + 2®K about 
15°C - f 13 „ c - 87.8 MHz. 

-Static frequency-temperature characteristic. 
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CONCLUSION 


The value of the dynamic T.C. a, obtained from the experimental frequen- 

_5 

cy responr.es is equal to 2.10 s/°K for ST-cut quartz SAW oscillator. This cof 
ficient is comparable to the "a of AT-cut quartz BAW oscillator (1.3.10 s/°K) 

But this dynamic thermal behavior can be modified by changing the mounting of 
the substrate in order to minimize the stress effects. 

The theoretical model must be improved : 

- by introducing the heat exchange by radiation 

- by using a bidimensional model which will take into account the edge 
conditions. 

This first study enabled to determine the order of magnitude of the SAW 
dynamic thermal effect and to make comparison with bulk waves. 

The continuation will consist in improving the theoretical model to study 
the influence of crystal orientation and wave propagation direction on the a 
coefficient and determine if doubly rotated cuts can be used with minimized dy 
namic effect. 
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